Spin rotation for ballistic electron transmission induced by spin-orbit interaction 
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We study spin dependent electron transmission through one- and two-dimensional curved waveg- 
uides and quantum dots with account of spin-orbit interaction. We prove that for a transmission 
through arbitrary structure there is no spin polarization provided that electron transmits in isolated 
energy subband and only two leads are attached to the structure. In particular there is no spin 
polarization in the one-dimensional wire for which spin dependent solution is found analytically. 
The solution demonstrates spin evolution as dependent on a length of wire. Numerical solution for 
transmission of electrons through the two-dimensional curved waveguides coincides with the solu- 
tion for the one-dimensional wire if the energy of electron is within the first energy subband. In the 
vicinity of edges of the energy subbands there are sharp anomalies of spin flipping. 
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PACS numbers: 72.10.-d, 72.25.-b 



I. INTRODUCTION 



The electron spin precession phenomena at zero mag- 
netic field induced by a variable spin-orbit interaction 
(SOI) in 2DEG systems was firstly proposed by Datta 
and Das [jjj as a way for the realization of the spin tran- 
sistor. For this, the spin precession is controlled via the 
Razhba SOI associated with the interface electric field 
present in the GaAs heterostructures that contains the 
2DEG channel g 



Vg = ha[p x a y -p y a x }. 



(1) 



The reason of spin precession is that the spin operators 
do not commutate with the SOI operator what leads to 
the spin evolution for electron transport. In particular 
the SOI has a polarization effect on particle scattering 
processes || and this effect was considered for different 
geometries of confinement of the 2DEG g |, g, @, g, §. 

The most simple case of the stripe geometry with the 
x-axis along the stripe and the z-axis perpendicular to 
the stripe gives the following transformation of spin state 
after transmission 



where § 



cos 0/2 
sin 0/2 



= 2m* aL 



(2) 



(3) 



and L is a length of the stripe. Therefore, the Razhba 
SOI induces a spin precession of the transmitted elec- 
trons. Notice that the spin precession is energy inde- 
pendent. This result is valid if the confinement energy 
h 2 /2m*d 2 , where d is a width of the stripe, is much larger 
than the spin-splitting energy induced by the SOI, and 
therefore, the intersubband mixing is negligible B. For 
the strong SOI the spin rotation angle becomes to be 
the Fermi energy depend for ballistic transport of elec- 
trons in the quasi one-dimensional wires and stripes g, ^) . 



Whatever the Razhba SOI leads to the spin precession in 
the (x, z) plane. Here we consider similar phenomena for 
electron transmission through the curved waveguide and 
quantum dots. The main difference between the straight 
waveguide and curved one is that the spin rotation is 
given by two angles. 

Next, we find out conditions under which there is no 
a spin polarization of transmitted electrons. We imply 
that a flow of incident electrons have no spin polariza- 
tion. By the spin polarization we consider the mean spin 
< a a >,a = x,y,z averaged over the electron flow. In 
particular for the transmission through quantum dot we 
show a principal role of the third lead for the spin polar- 
ization. 



II. THE SPIN-ORBIT INTERACTION IN THE 
INHOMOGENEOUS TWO-DIMENSIONAL CASE 

We write the total Hamiltonian of a confined 2DEG as 
h 2 ( d d 



H 



2m* 



d 

dx 2 



dy 2 



+ V(x,y) + V so (4) 



where V(x, y) is the lateral confining potential. Following 
to Moroz and Barnes [^| we assume that the SOI operator 
Vso is formed by three contributions 



Vso = VSl 



so 



V'ftQ 
so ■ 



The first Vg is related to the Razhba SOI ([!]), in which 
the SOI constant a proportional to the macroscopic 
interface-induced electric field is considered as constant. 
The second contribution Vg Q to the SOI comes from the 
electric field E(x, y) related to the confining potential. 

In order to derive the second contribution to the SOI 
we begin with general description of SOI |0 



Vso 



Am 2 c 



<7(Exp) + ^a(VxE)|. (5) 
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For microscopic electric field E the second term in (|5|) 
equals zero. However for model cases of the confining 
potential V(x, y) the electric field can violate an equality 
VxE = 0. In this case the second term in (||) is necessary 
to provide a hermiticity of the total SOI operator. 

For a 2DEG confined at semiconductor heterostructure 
interface we can reduce the z-coordinate performing aver- 
age over electron wave function tpo(z) strongly localized 
along the z-direction 

Vso => J dzMzWsoMz)- (6) 
As a result we obtain 

VgQ = — 7< <J Z {E X py - EyP X ) - E Z ((T x Py - <7yP X ) 



the SOI. Then we can specify the electron state by quan- 
tum numbers, the number of the energy subband n and 
spin projection a = a z . This assumption implies that 
far from waveguides or quantum dots the SOI constant 
a equals zero in the electrodes. Neglecting by real space 
behavior of the microscopic electric field at the edge of 
the heterestructure we assume that the field is directed 
normal to the plane of the heterostructure everywhere 
and has a stepwise behavior at the edges. As a result 
we obtain the stepwise behavior for the Razhba SOI con- 
stant a. Such a model was used by Hu and Matsuyama 
0. Similar to (@) we obtain that the third contribution 
to the SOI takes the following form 



—fkr z 
2 



' dE x _ dEy 
dy dx 



i u f dE z 



dE z 

T x 

dy 



(7) 

Here electric field components are considering in a mean- 
ing of integral (^) and depend on x, y only. 

For particular case of straight wire directed along the 
y-axis with the lateral confining potential U — U (x) 
we obtain from (R) the expression given by Moroz and 
Barnes (formula (5) in ||).They used a parabolic approx- 
imation for the confining potential. Here we consider a 
popular hard wall approximation and imply the following 
confining potential 



U(x) 



_/0, if \x\<d/2, 
U , if \x\ > d/2 



Then substituting the electric field E x = —U'(x) into (0) 
we have 



Vg (x) — hkja z sign(x)Uo5(x =F d/2). 



(8) 



For 1 2 1 > d/2 we have from the Schrodinger equation the 
following solution 



(9) 



where C is the normalization constant. Using a property 
of delta function that a difference between derivatives of 
the wave function at the right and left of the delta func- 
tion is obeying to A-0'(±d/2)) = ±2m*ka z -fU 'ilj{±d/2) 
we have from (|9|) that 

A^'(±d/2)) -> 

for Uq — > oo. Therefore in the hard wall approximation 
an effect of the second contribution Vg limits to zero. 

Next, for numerical computation of the transmission 
through the semiconductor heterostructure we assume a 
connection at least to two electrodes in which there is no 



III. THE TRANSMISSION THROUGH 
BILLIARD WITH THE SOI 



In this section we prove that the SOI gives no spin 
polarization for electron transmission through arbitrary 
billiards if energy of incident electron belongs to the first 
energy subband. In dimensionless form the stationary 
Schrodinger equation has the following form 

*-($;;])■ (id 



Here e = E/E ,E 



2m*L 2 ' 



VSO 



d_ 

dy 



i / d/3 
- I a-y 



dx 



dy 



(12) 

L is a characteristic scale of the system, (3 = 2m* aL is 
the dimensionless SOI constant. Correspondingly in Eqs 
( |TT| ) and (|T^) coordinates x, y are also dimensionless. 

Let S be an area of the structure under consideration 
which involves a billiard the SOI and leads as shown in 
Fig. 

Let r denote a boundary which crosses input lead and 
output one at Ti and T2 respectively. We suppose that 
there is no spin-orbit interaction in the leads, i.e (3 — 
at T{, i = 1, 2. At the rest of boundary T we imply the 
Dirichlet boundary conditions for solution of Schrodinger 
equation ( p]) i/j\ r =0. As the scale L we take L = d. 

Therefore we can write solution in the electrodes as 
follows 

\inc 7 n,a >= v2 sm(irny) exp(ik n x)\a >, 
\refl, n, a >= \[2 2J r mn ^^ sin(7rmy) exp(— ik m x)\a' >, 



\tr,n,a >= V2 2J tmn,ay smfrmy') exp{ik m x')\a' > 



(13) 



3 




FIG. 1: Schematical view of two-dimensional billiard with two 
attached leads. Dashed area shows a region with the SOI. The 
area S has boundary V which crosses input lead and output 
one at Fi and F2 respectively. 



where | a > is the spin states defined spin projection along 
some axis, say, the z-axis. The energy is 



/ 2 i 2 2 

e = fc„ + 7r n 



(14) 



where n — 1,2,... numerates a number of the energy 
subbands. 

Introduce complex derivatives 



d 1 / d .d 



dz 2 \dx dy / 
we write the Schrodinger equation (O) as follows 



(15) 



d 2 1 \ Pdu* ^W = () 
Ozdz* 47 Ul 2 dz 4 dz ' 

5^ + 4T 2 -2^"7a? = ' (16) 



where u\,U2 are the components of the spin state. As- 
sume that there is auxiliary degenerated state with com- 
ponents Vi,V2- In particular, it might be the Kramers 
degenerated state. Then, for these two states the Green 
formula follows 



J (uAv - vAu)dS = <j> 



r \ on on J 



where n is an exterior normal to the boundary F. 
From the Schrodinger equation we have 

1 \ /3 du 2 1 80 
dzdz* 4 / 2 dz 4 oz 

I \ P dv 2 1 9/3 n 
Ozoz* 4 / 2 az 4 



"1 




/ a 2 

a 2 1 \ p d Ul 1 a/3 

7 e U2 - o^iU-r " T^TTII = 0- (18) 



n p dvi 1 a/3 
Vazaz* 1 i e r 2 ~ 2 Ul a^~4 MlUl a^ = ' 



2 az* 



az* 



Combining each couple of equations in (18) we obtain 



''2 



Ul 



d 2 u\ 
dzdz* 

d 2 V 2 

dzdz* 



- u 2 



■vi 



d 2 v 1 



tf +-e(uiv 2 + U2V 1 ) + - 
dzdz* 4 2 



1 d(Pu 2 v 2 ) 



dz 



d 2 u 2 



dzdz* + r iuiV2+U2Vi) -2 



1 d{(3u\vi) 



dz* 



0, 



= 0. 
(19) 



Extracting the second equation from the first one in (J1J 
we obtain 



"2 



Ul- 



d 2 v 2 



U 2 - 



d 2 vi 



d 2 Ul 

dzdz* 1 dzdz* ' * dzdz' 
1 8(/3uiVi) 1 d(Pu2V 2 ) 



d 2 u 2 
dzdz* 



dz 



dz 



= 0. 



(20) 



Integration of this equation over the billiard area S 
with use of the Green formula dl^) gives the following 



dui 
r v an 



Ml—— I <M 

an 



a«i a^2 

r v dn dn 



d(2/3uivi) 
dz* 



dS 



d(2(iu 2 V2) 
dz 



dS = 0. 



dl 



(21) 



Since at T either ui = 0, Vi — , or /3 = 0, the last 
two integrals in ( pl[ ) equal zero and formula (|2l|) can be 
rewritten as follows 



E 

i=l,2 



E 

»=1,2 



I'2" 



a?xi 



Ul 



dv 2 
dn 



«2 



a^i 

dn 



Vl 



du 2 
dn 



dl 



dl = 0. 



(22) 



This formula is sufficient to establish some symmetry 
rules between ingoing and outgoing states. Let us con- 
sider the first-channel transmission for e < 47r 2 . In order 
to ignore evanescent modes we will consider that bound- 
aries Vi cross the leads far from the scattering region as 
shown in Fig. [I] . Let electron incidents from the in- 
put lead being completely spin polarized up. It means 

/ 

note the corresponding state interior the structure S as 



that for the incident state 



a >-- 



We de- 



which is used as the u-solution in Eq. (|22j). 



Correspondingly 



, denotes the v-solution in 

U2i(x,y) _ 

Eq. ( p2| ) for the case of electron incidenting with spin po- 
larized down . We suppose that the boundaries Ti and 
T2 cross the leads normally the leads and the x-axis is 
parallel to the leads. Hence the normal n is parallel to 
the x-axis. Then from ( |l3] ) at the boundary T2 which 
crosses the output lead we obtain the following relations 



df 
dn 



ikif 



(23) 
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where function / refers to all components 

These relations allow to exclude the boundary r 2 from 
(^||). At the boundary Ti which crosses the input lead 
we have 



dn 



i/cittij — 2ik\ sin(7ry), 



du 



21 



dn 

dun 

dn 



= ikiu 2 -\, 



du 



21 



On 



ikiU2i — 2ik\ sin(7ry). 



(24) 



We imply here that the origin of the x, y coordinate sys- 
tem is at the boundary T\. Substituting the relations 
(|J) into Eq (H) we obtain 



(wi T - u 2 i) sin(7ry) dy = 0. 



(25) 



Since at the boundary T± 

u n = ui T (a;)sin(7ry),'U2i = u 2L (x) smfry). 
we obtain from ( p5| ) 

"it = u 2i- 

Thus from ( p"3| ) it follows that amplitudes of the reflection 
r T,t= r W- ( 26 ) 
Next, we take that the state ' Ul ' — ' 



cides with the state 
simplifies as follows 



E 



U 2 J \U 2 -\ 

in (H). Then equation (|f) 



u 2 



dn 



dn 



i=l,2 Fi 

Substituting into this formula relations (^ij) we obtain 



2ik I U2] sin(7ry) dy = 0. 
'r. 



(27) 



It gives us that U2] = or according to ([13]) r-fj = . 
Also similarly we obtain that = at the boundary 
T\. Thus we can write the second symmetry rule for 
reflection amplitudes 



"t,i = n,t = o. 



(28) 



From symmetry rules ( |26| ) and (28) and from the current 
preservation it follows that the transmission probabilities 



(29) 



do not depend on the spin polarization of incident elec- 
tron. 

Till now we considered incident waves as spin polarized 
along the z-axis at the boundary Ti. Let now consider 
a flow of incident electrons which have no averaged spin 
polarization. In particular we can present that half of 
electrons have the incident state with spin up and half of 
electrons have the incident state with spin down. Let us 
consider corresponding transmitted waves at the bound- 
ary r 2 . We prove that for a transmission through the 
billiard with two attached leads there is no averaged spin 
polarization, i.e. < a a >= 0, a = x,y, z if electron inci- 
dents being spin unpolarized in the first energy subband. 
As previously we take the incident state in the form (|lF 
and write the states in leads as 



|V T >= 



wit 

u 2] 



Hi >- 



U21 



(30) 



where the arrows up and down indicate that electron in- 
cident with spins up and down. We take in the Green 
formulas ( p2| ) the first function u as > and the sec- 
ond function v as a y C\ipi > where C means a complex 
conjugation. It means that the second function is the 
Kramers degenerated state. Hence 



_ I «1T 

W2T 



2 i 

-iu" 



1! 



(31) 



Let us calculate integral (J22j). At a boundary Fi cross- 
ing the input lead the integral equals zero since in the 
input lead U2 = 0, v\ = 0. The second contribution into 
integral (^2|) relates to the boundary F2 crossing the out- 
put lead. Using transmitted solution ( ft3| ) one can write 
Eq. (|H) as follows 



dux dv 2 . 

u i-^~ 1 dl 

dn dn 



dvx du 2 ■ ,, 

"2^ v i^r~ I dl 

r \ dn dn 



= 2ik d> (u 2 ui — U2V\)dy = 0. 

Therefore U1V2 = u 2 «i, or in terms of notations ( |3l"| ) 

wiT"ix = - u 2t u 2i- (32) 
From ( |32j ) it obviously follows 

KtlKJ = l«2tlKil- 

Moreover relation ( p9| ) implies that 

| UlT | 2 + | M2T | 2 = | Uli | 2 + | U2i | 2 
From these two relations one can obtain that 

I«2tI = |ui;|,Ktl = KJ- (33) 
Finally relations (H|) and (|3|) give 

"iT u 2T = - u n u 2i- (34) 
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FIG. 2: Spin polarization of electrons transmitted through 
the three terminal quantum dot versus energy of electron in 
the first energy subband. An inset above shows a geometry 
of the structure. 
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Mean values of spin components in corresponding 
states ( |30| ) are the following 

< <J X >t= i?e(uifu| T ), < a y > T = 7m(Mi T U2 T ), 

< (T 2 >|= |ui T | 2 - |u 2 t| 2 , 

< °x >x= Re{u n u* 2l ), < a y >;= /mfMijiiy), 

< er z > x = |u u | 2 - |u 2i |- (35) 



(36) 



Eqs (Bah - (Ba) give rise to 



< cr Q > T = - < a a >i, a = x,y,z, 



i.e. the spin polarizations are exactly opposite in sign 
for transmission of electrons incidenting in corresponding 
spin polarized states. 

Thus, for the transmission through any billiard with 
the SOI with two attached leads the spin polarization 
does not exist if the flow of electrons incidents in the first 
energy subband and have no spin polarization. Also, if 
there is no intersubband transmissions t mn . aa i =0,m/ 
n, the spin polarization equals zero for arbitrary energy. 
It takes place approximately, for example, for adiabatic 
structures similar to curved waveguides (section V) . How- 
ever in a vicinity of edges of the energy subbands 7r 2 n 2 
the SOI gives rise to intersubband mixing. As a result we 
obtain in numerical calculations strong spin polarization 
near the edges. Moreover, if the billiard is connected to 
three or more leads, the spin polarization of transmit- 
ted electrons exists even for the transmission in the first 
energy subband. The effect of the third lead is demon- 
strated in Fig. ^. Hence this effect propose a way of 
the spin transistor complimentary to the way proposed 
by Datta and Das jjj . The spin polarization of transmit- 
ted electrons can be governed by a value of connection 
of the third lead with the quantum dot. The most sim- 
ple way is to apply local electric field in the vicinity of 
the connection which implies potential barrier closing the 
connection of the dot with the third lead. 



FIG. 3: Schematical view of one-dimensional curved wire. 

IV. THE ONE DIMENSIONAL CURVED WIRE 

A model in which only the single channel transmission 
takes place is the one-dimensional wire. Therefore for 
a transmission through the one-dimensional wire of any 
form the SOI can not give rise to the spin polarization. 
However this model is interesting by that allows to find 
spin evolution analytically. A case of straight wire was 
considered by [Q, [|. Here we consider a curved wire 
consisted of a segment of circle with radius R attached 
to infinite straight one dimensional wires as shown in Fig. 

1 

We take a length of the segment as L = 4>qR and a 
position coordinate as s = 4>R. The Hamiltonian of the 
wire has the following form jl2|, [l3j 



H = 



2m* R 2 



H. 



H = 



<j T sin ( 



2 

T< (37) 



where (3 — 2m*aR is the dimensionless SOI constant. 
Since [J Z ,H] = where J z = ~^~§^ ~ ^z, a particular 

solution of the stationary Shcrodingcr equation H\if) >= 
e\ip > has the following form [jl^, [l3|, |l4| 



\4>>= 



(38) 



The parameter /i defines the dimensionless wave number 
as k = (i/R and is arbitrary until the boundary condi- 
tions are imposed. Substituting the state ( |3^ ) into the 
Shcrodinger equation one can obtain the following rela- 
tion between the energy of electron e and the wave num- 
ber jj, 



(e- )U 2 )(e-( M -l) 2 )- y 5 2 ( / i-l/2) s 







(39) 
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The matrix U{<j>) has the following form 





U{<j>) = f ex PW/ 2 ) 



exp(— i(j>/2) 



(45) 



Evolution of the electron state (|42|) as length s = <pR 
of the curved wire can be presented as 



|^(0) >= e^+^-^/^^A 




FIG. 4: The energy spectrum defined by formula ( fio| ) for 
a — 1 . Values of /i corresponded to clockwise movement of 
electron along the curved wire are shown by thick points. 



which gives 



e„ = (/i - 1/2) 2 + 1/4 + v\ii - \l1\yffP + \Ji, v = ±1. 

_ _ (40) 

The spectrum ( (40[ ) is shown in Fig. Q. 

For fixed energy e Eq.([i"o|) gives four solutions for the 
wave number /i. It is well known that || for electron 
transmission through potential profile a reflection is neg- 
ligibly small if the characteristic length of inhomogeneity 
much exceeds the wave length (adiabatic regime). For 
our case we assume that the radius of curvature of the 
wire is much larger in comparison with the electron wave 
length. So we can ignore the reflection for electron trans- 
mission through the quasi one-dimensional waveguide. 

Since there is no reflection for transmission through 
the one-dimensional waveguide we need only those val- 
ues of the wave number /i which correspond to clockwise 
movement of electron in the waveguide. We denote its 
as /ii,/i2 shown in Fig. |J . In what follows we use the 
following relation between \i\ and fi 2 : 



2A 



A*2 - Mi 



V1 + /3 2 



(41) 



Therefore general solution of the Shcrodinger equation 
for the electron transmission without reflection can be 
written as follows 



m<t>) >-- 



v=\.1 



U(<t>) 



where 



-i/3 



i(3 -1 -y/l + 2 



-i/3 



1 + ^T+W 

i!3 -1 + ^TT]P 



A 2 
B 2 

A x 
Bx 



(42) 



0, (43) 



= 0. (44) 



exp(iX<j>) 






exp(— 



where 



A 



A 2 A x 
B 2 Bx 



A"V(0) > (46) 



(47) 



From Eqs (|43|) and ( J44| ) we can rewrite (|47|) as follows 

A 2 -BX 



A = 



and 



B 2 A 2 



i(3 



Bo 



(48) 



(49) 



Eq. @ can be presented as \ip(cj>) >= T(<p)\i/)(0) > 
which shows that the unitary matrix T((f>) has a meaning 
of the transfer one. Since the state \ip((f>) > is spinor 
one the matrix T((f>) corresponds to rotation matrix for 
transport of electron along the wire. In general case the 
rotation matrix is given by the Euler angles {ip, 9, 7) |lEj ] 

R{<pAl) = e- l ^°>e- l6 ^ a *e- l ^ a >. (50) 
The rotation matrix R has the following form 

_( exp(— i\(tp + 7)) cos(^6*) — exp(— \i{<p — 7)) sin(i^ 
~ ^ exp(ji(v?-7))sin(i6>) exp(i±(<p + 7)) cos(±0) 

(51) 

In order to find the Euler angles let us consider to which 
rotation corresponds matrix 



(52) 



exp(iA0) 
exp(— iXcf)) 

If the matrix A were unit, the matrix 



exp(iA0) 
exp(— i\(j>) 



would correspond to the rotation by the angle 9 = —2X<fi 
around the z-axis. The matrix A in J52I ) gives rise to the 
clockwise rotation around the x-axis by the angle 9 which 
satisfies to the following equation 



^ = cot(0/2)e-^ = % 



(53) 
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FIG. 5: A fragment of two-dimensional curved wire with 
width d — 1. 



where 

TJ A -1/2 ^ -1/2 <5 -1/2 <9 1/2 <9 [rr\ 

and d is a width of the waveguide. In what follows we 
consider a segment of the two-dimensional ring with con- 
stant curvature 7 = 1/R attached to straight leads with 
the same width as shown in Fig. ||. Therefore for the 
segment we can write 

x = a(s) — ub'(s) 
y = b(s) + ua'(s) 
a(s) = -Rcos(s/R),b{s) = Rsin(s/R) 

g ^ = l + ul { S ) = ^- (56) 

with 7(s) as a curvature of the curved waveguide which is 
taken below constant. The SOI takes the following form 
at the curved part of the waveguide 

VSL = 



In order to fulfill this equation we choose 



(f = 7r/2, cot • 







(54) 



The angle 7 is remaining undefined. Below we put 7 = 0. 
Let us choose new axis z' in the (y, z) plane with the angle 
as the angle between the z-axis and z'-axis. Thus, the 
full rotation matrix consists of the anticl ockwise rotation 
around the the z'-axis by the angle (f>\/l + (3 2 and the 
clockwise rotation by the angle <fi around the z-axis. The 
last statement follows from the matrix U ((f)) in Ec 
A knowledge of the evolution of the spin state 
dependent on the length <j)R of the curved waveguide 
gives us a possibility to calculate in particular evolu- 
tion of the spin components. The result of calculation 
is shown in Fig. ^ by squares, triangles and circles for 
spin components u z , a x , a v respectively. As seen the re- 
sults of numerical computation for the two-dimensional 
curved waveguide are surprisingly close to present one- 
dimensional model. 




^(lr + *2- 1/2 £) 







(57) 

At the leads we assume that there is no the spin-orbital 
interaction (/? = 0) as well as 7 = 0, g 1 ^ 2 = 1. 
The Shcrodinger equation 



with the total Hamiltonian as H = Ho + vso takes the 
following form 

1/2 ± f.-Vi^ +fl -i/al ( n ^^l\ 



ds \ ds J ^ du \ du J 



1/2 !L („-v2f£tk\+„-v*JL („v**tk 



ds \ ds ) ® du \ du 



V. THE TWO-DIMENSIONAL CURVED 
WAVEGUIDE 

For consideration of the two-dimensional curved 
waveguide we introduce the curved coordinate system 
(s, u) jl6[ [lTj where s is the coordinate of central line 
along of the waveguide shown in Fig. ||. 

We express the Hamiltonian of the waveguide in di- 
mensionless form by following way 



+ «* + *^«(£-*-"g)-0, (5S) 

The solutions of Eqs ( j58[ ) which satisfy to the Dirichlet 
boundary conditions (u — ±1/2) can be presented as 




V>t(w, s) = ^ A*\ n (s) sm(Trn(u + 1/2)) 

n=l 
00 

^(u, s) = A ln {s) sin(7rn( U + 1/2)). (59) 

n=l 
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Substitution of fl59| ) into Eqs (|58|) gives 

oo 

^2[L mn A" n (s) + P mn A u {s) - (3e^ s Q mn A Ln (s) 
-iPe^ s R mn A' ln (s)} = [(Trm) 2 - e]A ]m 

OO 

n=l 

- z/3e-^ s i? m „A' T „( S )] = [(^m) 2 - e]^ im . (60) 
Here we introduced the following notations 

f 1/2 sin(7rro(u + 1/2)) sin(7m(u + 1/2)) 

-^mn — * / 7^ — ; T5 MU) 

J-i/2 (I + U7) 2 

/" 1/2 7sin(7rm(M + l/2))cos(7rn(M + 1/2)) , 
r mn = 27m / — du, 



Rm.n. 2 



1/2 1 + «7 

1/2 sin(7rm(u + l/2))sin(7m(u + l/2)) 



-1/2 
1/2 



1 + 1*7 



27m / sin(7rm(u + 1/2)) cos(7rn(u + l/2))du. 

-1/2 



-t — 1 

CD uo 
C 

o 

CL 

O 
O 

"I.-0.5 
CO 



^vavs4 v e 



® V 



(61) 




VI. NUMERICAL RESULTS 

In numerical practice we solve the system of Eqs ( |60| ) 
and ( |6l"l ) taking a finite number of waveguide modes. This 
number of modes was controlled by the normalization 
condition that sum of the total reflection probabilities 
and the total transmission ones is to be equaled to unit. 
The spin components < cr a > were calculated at the 
attached outgoing straight electrode in which we assumed 
there is no the spin-orbit interaction by following formula 



1/2 

J_' /2 du < ip{u, s)\ap\ip(u.s) > 
< apis) >= !— . 

I-l/2^ U ^ 1p{v-i s)\lp(u.s) > 



(62) 



In Fig.|B| the outgoing electrode as well as incoming one 
are not shown. Fig. ra shows evolution of the spin com- 
ponents ( p2| ) versus the longitudinal coordinate s. 

It is surprising that for energy of incident electron 
far from the edge of energy subband the spin evolution 
almost coincides with the one-dimensional curved wire 
shown in Fig. |^ by squares, triangles and circles. In 
Fig. @ (a) the energy dependence of the spin components 
are shown which demonstrates remarkable phenomenon 
of spin flipping at the edge of the second energy subband 
E-x = (27r) 2 39.4. It is interesting that increasing of re- 
gion with the SOI by increasing of length curved waveg- 
uide or increasing of the spin-orbit constant leads to dou- 
ble flipping of electron spin for transmission through the 



FIG. 6: The spin components as dependent on the length s 
of a curved two-dimensional waveguide. The result of calcu- 
lation based on the state (^) for the curved one- dimensional 
wire is shown by squares (er z ), triangles (a x ) and circles (a y )- 
The radius of wire R = d where d is the width of waveguide. 
The dimensionless spin-orbit constant /3 — 2m* ad equals 
unit. (a) The dimensionless energy e = 25 (the first channel 
transmission) and (b) e = 39.25 (near an edge of the second 
subband). 



waveguide as shown in Fig. |?] (b) . This phenomenon is 
a consequence of the intersubband mixing by the SOI as 
it was discussed in section III. 

Therefore one can expect strong deviation of the 
curved two-dimensional waveguides from the onc- 
dimensional one for the spin evolution near edges of the 
subbands 7r 2 rz 2 . In fact, one can see from Fig. |7| (b) that 
for energy of incident electron E w 47r 2 the spin evolution 
with length of the curved waveguide is strongly deviates 
from the case of one-dimensional curved wire. 



Acknowledgments 

This work has been partially by RFBR Grant 01-02- 
16077 and the Royal Swedish Academy of Sciences. 
* e-mails:almsa@iffn. liu.se, almas@tnp.krasn.ru 



9 





FIG. 7: The spin components as dependent on the energy of 
incident electron for (a) <j>o — 90°-curved waveguide and (b) 
cj>o = 180° one. = 1. 
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